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ABSTRACT 

We use cosmological SPH simulations to study the cool, accreted gas in two Milky Way-size galaxies through cosmic 
time to z = 0. We find that gas from mergers and cold flow accretion results in significant amounts of cool gas in 
galaxy halos. This cool circum-galactic component drops precipitously once the galaxies cross the critical mass to form 
stable shocks, M v ; r = M s h ~ lO 12 Af . Before reaching M s h, the galaxies experience cold mode accretion (T < 10 5 
K) and show moderately high covering fractions in accreted gas: f c ~ 30 — 50% for R < 50 co-moving kpc and 
Njjj > 10 16 cm~ 2 . These values are considerably lower than observed covering fractions, suggesting that outflowing 
gas (not included here) is important in simulating galaxies with realistic gaseous halos. Within ~ 500 Myr of crossing 
the M s h threshold, each galaxy transitions to hot mode gas accretion, and f c drops to ~ 5%. The sharp transition in 
covering fraction is primarily a function of halo mass, not redshift. This signature should be detectable in absorption 
system studies that target galaxies of varying host mass, and may provide a direct observational tracer of the transition 
from cold flow accretion to hot mode accretion in galaxies. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: evolution — galaxies: halos 
— methods: numerical 



1. INTRODUCTION 

Recent advances in galaxy formation theory have em- 
phasized the importance of cool gas accretion onto 
galaxies: gas that never shoc k-heats to the viria l tem- 
perature of the halo (e.g.. IBirnboim fc Dekell 120031 : 
IDekel fc Birnboiml l2006t IKeres et all 120091) . Under this 
picture of galaxy formation, the cooling time of gas enter- 
ing halos less massive than a critical threshold, M v i r = 
M s h, is too short to sustain compressive shocks, with 
M 8 h ~ 10 n - 5 ~ 12 M Q . Galaxies within dark matter halos 
less massive than M s h experience "cold mode" accretion, 
since baryonic accretion onto these galaxies is dominated 
by cool gas. Galaxies above this transition experience 
"hot mode" accretion, with infalling gas shock-heating 
to the virial temperature. In practice, the precise value 
of M s h depends somewhat on definition. In what follows 
we are interested in the absolute shut-down mass, when 
almost no cool accreted gas reaches the central regions 
of the halo. Thus, we adop t M s h ~ 10 12 M Q , motivated 
by IDekel fc Birnboiml i{2()0(i: for shocks at large fractions 
of the virial radius for Z ~ 0.1 metallicity gafl 

Unfortunately, there are currently no definitive ob- 
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servational tests to detect cosmological cool gas accre- 
tion. On the contrary, numerous observational studies of 
cool halo gas around galaxies have e mphasized the pres 
ence of gas ou t flows, not inflows (e. g . | Steidel et al.lll996l 
Shaplev et all 120031: iWeiner et all 120091: iSteidel et al 
20101: iRubin et al.l \20lW . The stark contrast between 



theory and observations is understandable at high red- 
shift, as gas inflow to galaxies at z > 2 is expected to flow 
along den se filaments, resulting in small g l obal coverin; 
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fracti ons (jFaucher-Giguere fc Keresl 12011 
12010(1 . In addition, galaxies at z ~ 2 are at the peak of 
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cosmological star formation (jHopkins fc Beacoml l2006h : 
one might expect feedback processes to dominate any 
observational indicator of gas accretion at these epochs. 

At lower redshifts star formation rates decline and the 
gaseous halos of galaxies are observe d as quasar absorp- 
tion s y stems (e.g..lBergeron fc Bo issc 19 91t iBowen et all 

' Barton fc Cookd 120091: 
These stud- 
ies are primarily of metal lines (Mgll, CIV, OVI, etc.) 
though with COS on HST we expect Lyman a observa- 
tions to increase. In this letter, we utilize cosmological 
hydrodynamic simulations to study cool gas accretion 
and the possibility of detecting it as quasar absorption 
systems. To correctly produce metal lines in a simulation 
requires radiative transfer, metal diffusion and modeling 
of local ionizing sources; however, we are focusing on the 
qualitative behavior of halo gas, so we will instead give 
results in terms of HI column density calculated in the 
optically thin limit, without local sources. For column 
densities below the Lyman limit (2 x 10 17 cm -2 ) this 
should be fairly robust, but for higher column densities 
we expect a full treatment would lead to quantitative but 
not qualitative differences. 

2. SIMULATIONS AND ANALYSIS 

2.1. The Simulations 

Our two simulations utilize separate sets of cosmolog- 
ical initial conditions, with each simulation tracking the 
evolution of a roughly Milky Way size dark matter ha- 
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los until z = 0. We refer to these two simulations as 
"Halo 1" and "Halo 2", since we primarily investigate 
the properties of the single most massive galaxy in each 
simulation. Halo 1 has an active merger history until 
z ~ 1.5, but subsequently experiences a relatively quies- 
cent merger history (a dark matter only simulation of the 
same initial conditions was performed at very high reso- 
lution in the Via Lactea II simulation of iDiemand et aT] 
(|2008| ) In contrast, Halo 2 experiences a more quiescent 
early history, but has a major merger at z ~ 0.8. For 
both simulations, the most massive galaxy has a dark 
halo mass of M v j r (z = 0) = 1.4 x 10 12 Af Q , and is posi- 
tioned within a large scale filament of dark matter and 
gas, typical for galaxies of this mass. 

We use multiple mass particle grid initi al conditions 
genera ted with the GRAFIC-2 package ([Bertschingerl 
(|2001| )) and the best-fit cosmologi cal parameters of th e 
WMAP three-year data release (|Spergel et all l2007h : 
Om = 0.238, fl A = 0.762, H = 73km s" 1 Mpc" 1 , n s = 
0.951, and erg = 0.74. We implement the "zoom in" tech- 
nique to prope rly account for large- scale tidal torques 
as described in iKatz fc White! (|1993f ), and each simula- 
tion is contained within a periodic box of 40 co-moving 
Mpc, with the highest resolution region limited to a ~ 6 
co-moving Mpc cube. We use the sm oothed particle 
hydro dynamics (SPH) code GASOLINE (jWadslevet. all 
120041 ). In the highest resolution region of each simula- 
tion, the masses of the simulated particles in the initial 
conditions are: (mdark, rn gas ) = (17,3.7) x 10 Mq, with 
a force softening of 332 pc. 

The code assumes a unifor m UV backgroun d from 
QSO, implemented following lHaardt fc Madaul (|1996f ) 
and F. Haardt (2002, private communication). It imple- 
ments star formation, as well as Compton and radiative 
cooling, as described in IKatz et ail (jl996t ). calculating 
the abundance of neutral hydrogen by assuming an op- 
tically thin ideal gas of primordial composition and in 
ionization equilibrium with the UV-background, treat- 
ing collisional ionization, photoionization and recombi- 
nation processes. We do not include full treatment of 
radiative transfer. Because we primarily focus on quali- 
tative attributes of cool halo gas, and its evolution with 
time, we do not believe this will significantly impact our 
results. We expect that a full treatment would increase 
the amount of very high density gas, causing a progres- 
sively larger increase to our reported covering fractions, 
starting at a few times the Lyman limit. 

Type II supernovae are modeled using an analytical 
treatment of blastwaves, creating turbulent motions in 
nearby gas particles that ke eps them from coolin g and 
forming stars, as described in lStinson et. a l lpOOfih . This 
feedback model results in minimal winds of ~ 100 km/s 
that mostly affect hot gas and are mor e prominent in 
M vir < 10 n M Q halos (jShen et al.ll20Tol) . For the mass 
range considered in this letter, this feedback model does 
not result in cool outflows that would be detectable 
through absorption. The only two free parameters in the 
star formation and feedback model have been fixed in 
order to produce galaxies with a realistic star formation 
rates, disk thicknesses, gas turb ulence, and Schmidt law 
over a range in dynamic masses (|Governato et al.ll2007l ). 
We note our feedback model is very similar to those used 
in recent simulations that have shown great success in 



produ cing realistic disk-type galaxies dGovernato et al.l 
I2009D. as well as m atching the mass-metallicity relation 
( Brooks et al.ll2007l ) and t he abundance of Dam ped Ly- 
man a systems at z = 3 dPontzen et al.ll2008h . We re- 
fer the reader to iGovernato et al.l ()2009l ) for a more de- 
tailed description of the simulation code. At each output 
snapshot of our simulati on, we define the virial radius by 
iBrvan fc Normanl (|1998l ) , noting that this is a fairly typi- 
cal definition used in iV-body simulati ons of dark matter 
substructure (e.g. iStewart et alll2008| ). 

While a detailed morphological analysis of our sim- 
ulated galaxies is beyond the scope of this letter, visual 
inspection shows that both galaxies are disks before they 
reach the transition mass for hot mode accretion. Af- 
ter this transition, Halo 1 remains roughly disk-like but 
grows a massive bulge as hot gas cools from the halo, sim - 
ilar to the most massive g alaxy in lGovernato et all ()2007l ) 
and iBrooks et all (j2009f ). Due to the major merger in 
Halo 2 just prior to reaching the transition mass, it de- 
velops into a more spheroidal galaxy. We note that pre- 
vious works that have produced disk galaxies at z = 0, 
even for galaxies with major mergers at z < 1, typically 
involve halos less massive than those we study here (e.g., 
7 x 10 n M Q in IGovernato et al.l <p009[) ). which may be 
the cause of these morphological differences. 

2.2. Analysis 

To study the cool halo gas in our simulated galaxies, 
we analyze 794 regularly spaced sight lines within 100 
co-moving kpc, for three orthogonal orientations of each 
galaxy. For each sight line, we compute the column den- 
sity of neutral hydrogen per unit velocity, as a function of 
velocity along the line of sight, mimicking the practice of 
searching for absorption gas along a quasar line of sight. 
We do this by dividing each sight line into spatial bins of 
order the size of the spatial resolution of the simulation. 
For each spatial bin, each gas particle with a smoothing 
length that intersects the line of sight contributes to the 
total mass density of the bin. For each of these parti- 
cles, the contribution of the spline kernel to the column 
density of each bin that overlaps the kernel is integrated 
along the line of sight. The velocity of each bin is the 
mass-weighted velocity of particles contributing to the 
total column density of neutral hydrogen. 

We define the covering fraction of accreted neutral 
hydrogen, f c (< R), as the total fraction of sight lines 
(within a radius R from the center of the galaxy) for 
which N m > lO 16 " 18 - 7 cm" 2 . For > 95% of cases, the 
full width half maximum spread in velocities of absorbing 
sightlincf0 is < 110 km/s. We choose these column den- 
sity limits in order to probe different types of cool gaseous 
structures that may exist in the halos of galaxies, noting 
that they correspond to a plausible minimu m threshold 
for halo gas detection via M g II absorption (|Rigbv et al.l 
[20021: [Churchill et al.|[2000h . 

Some past studies of metal line absorption systems as- 
sociated with galaxies suggest that there may be a critical 
radius, outside of which cool gas clouds can not form (e.g. 
iChen et al.ll200ll [20To| ). while a study by iSteidel et all 
(|2010D suggests that f c (<R) decays as a power law in 

9 Because our feedback model docs not produce cool gas out- 
flows, our absorption sightlines are never offset from the systemic 
velocity by more than 1, 000 km/s. 
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Fig. 1. — Masses and covering fractions of each galaxy versus look-back time. The left panels show Halo 1, while the right panels show 
Halo 2. Top: the total virial mass (solid black), and baryonic masses within R < R vlT /3 (~ 100 kpc at z = for halos of this mass). Stellar 
mass, hot gas mass (T > 10 5 K), cool gas mass (T < 10 5 K) and neutral hydrogen are given by the dashed green, dash-dot red, dotted blue, 
and dotted magenta lines. The thick horizonal lines on the right side of each plot show the cool gas mass and neutral hydrogen contained 
within the central galaxy and massive satellite galaxies at z = 0. Bottom: covering fraction of neutral hydrogen within R < 50 co-moving 
kpc as a function of time. Black, red and blue lines represent minimum column densities ./Vjji > 10 16 0,17 - 3 '- 18 ' 7 cm -2 , respectively. The 
vertical hashed bar in each panel shows when the galaxy transitions from cold mode to hot mode accretion (see Q . 



R. Based on the typical extent of neutral hydrogen in 
our galaxies, we present covering fractions in terms of 
three fixed radii: R < 30, 50, 100 co-moving kpc (for a 
det ailed exploration of the radial dependence of f c (<R), 
see lStewart et al.ll2011[ ). We include the galaxy disk in 
our analysis, but its contribution to the covering fraction 
is relatively minor. For R < 30, 50, 100 co-moving kpc, 
f c for the disk is < 10, 4, 1%, respectively. 

3. GALAXY GROWTH AND EVOLUTION 



Figure Q] shows the mass accretion history of our galax- 
ies, as well as the covering fraction of each galaxy over 
time (Halo 1 on the left, Halo 2 on the right). In the top 
panels, the virial mass is given by the solid black line, 
while other curves give baryonic masses within i? v j r /3 
(100 kpc at z = 0). The stellar mass, hot gas mass 
(T > 10 5 K), cool gas mass (T < 10 5 K), and neutral 
hydrogen mass are given by the dashed green, dot-dashed 
red, dotted blue, and dotted magenta lines, respectively. 
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Halo 1 experiences several large mergers at z > 1.5, with 
a relatively quiescent growth after z = 1, while Halo 2 
has a more quiescent early growth (z < 1), but experi- 
ences a major merger at z ~ 0.8, and another moderately 
large merger at z ~ 0.2. These signatures, especially the 
late-time mergers for Halo 2, can be seen by the sharp 
increases in various curves in the top panels. Note that 
most of the halo gas within R < i? v ir/3 is cool until the 
halo grows to M v ; r <~ 1O 12 M . (See $4] for more on this 
transition. The vertical hashed bar in each panel shows 
the range where M v i r = 0.9— 1.0 x 10 12 Af© for each 
galaxy.) For proper context, we also include the cool gas 
mass enclosed within the main galaxy (R < 10 kpc) or 
in massive satellite galaxies at z = 0, represented by the 
thick horizontal lines at the right side of each plot. It 
is clear, then, that the vast majority of cool gas is as- 
sociated with galaxies at late times and is not spread 
throughout the halo. In contrast, when the galaxy first 
reaches the transition mass, the cool gas in the galactic 
disk is only ~ 10 9 Af Q , and does not dominate the cool 
gas mass within i? vu -/3. 

The bottom panels of Figure Q] show the time evolution 
of f c (R < 50 co- moving kpc) for each galaxy, with dif- 
ferent curves representing different column densitieJ^I: 
Njji > 10 16 cm" 2 (thick black), N HI > 10 173 cm" 2 
(thin red), and Nhi > 10 187 cm -2 (thin blue) In general, 
both galaxies show remarkably flat evolution of covering 
fraction with time (for M vir < A/ sh ), with f c ~ 30 — 40% 
for z < 4. In detail, each galaxy's covering fraction varies 
with recent merger history, as demonstrated by the sud- 
den spikes and dips in Figure [U These variations are pri- 
marily caused by mergers, which play an important role 
in cumulative gas accretion onto galaxies as well as halos 
(especially at z > 1 , when mergers are mor e frequent and 
more g as-rich than at late times; see e.g.. [Stewart et al.l 
I2009afbh . We note that IKacprzak et al.1 (|2010l ) have 
carried out a similar investigation to what we present 
here, and also found an accretion-related origin for many 
Mg II absorption systems. In addition, IKacprzak et al.1 
(|2007|) found a correlation between galaxy asymmetry 
and the presence of halo gas, again suggesting a link be- 
tween gas accretion and halo absorption systems. Re- 
assuringly, our covering fractions at high redshift are 
in qualitative agreement with a hig h- redshift simulation 
by iFaucher-Giguere fc Keresl (|2011| ) , utilizing a different 
SPH code and including treatment of radiative trans- 
fer. For similar mass galaxies as our own, they quote an 
average covering fraction of ~ 20 — 30% at z = 2 — 4 
(R < i?vh/2 and N m > 10 172 cm~ 2 ). 

We emphasize that these moderate covering fractions 
are a result of gas accretion only. Our feedback model 
does not produce significant outflows. We do this both 
to achieve a robust result (the physical mechanism be- 
hind galaxy outflows is largely uncertain) and to focus 
on signatures of cold flow accretion. Since galactic out- 
flows should populate the halo with additional cool gas, 
our results serve as a lower limit on the true covering 
fractions of galaxies. 



10 Note that the higher column density lines will not be correct, 
in detail, because we do not have full treatment of radiative transfer 
in our analysis. Our estimates are likely too low for higher density 
gas. 



4. COLD FLOWS AND COVERING FRACTIONS 

Ob served covering fractions typically vary from ~ 5 0% 
(e.g., iTripp fc Bowenl l2005t iChen fc Mulchaevl [200l t o 
near unity (e.g.. iChen et al. l l200lHChen fc Tinkerll200"cl . 
but there are indications that more massive galaxies 
(LRG s) may have systema t ically lower values (~ 10 — 
15%, IGauthier et ail l20Tot 1 owen fc Chelouchjl2l)T0h . 
Despite this variation in the literature, our findings are 
significantly lower than those observed, suggesting that 
outflows are a required component in producing galax- 
ies with realistic gaseous halos. Still, the contribution 
from accreted gas presented here is significant, so that 
the sudden drop in f c (Figure [1]) should be observable in 
real galaxies. 

Keeping in mind the theoretical framework of cold and 
hot mode accretion (see ?JTJ) , Figure [2] shows the covering 
fraction as a function of halo virial mass rather than 
time, for three choices of outer radii: 30, 50, and 100 
co-moving kpc (dot-dashed blue, solid black, dotted red 
lines, respectively). Regardless of the radius used, f c (< 
R) declines sharply at M v j r = M sri , when cold mode 
accretion ends for each galaxjEH- 

A preliminary analysis suggests the timescale for cool 
halo gas to fall onto the central galaxy is short (Stcw- 
artllc et al. in preparation); thus the reason for this 
sharp decline in covering fraction becomes clear. Once 
the galaxy crosses the threshold into hot mode accretion, 
the existing cool halo gas quickly sinks onto the galaxy 
and cannot be replaced by subsequent gas accretion. A 
galaxy in this hot mode accretion can only maintain a 
cool gaseous halo by mechanisms other than direct ac- 
cretion. For example, outflowing winds may be capa- 
ble of re-distributing cool gas from the galaxy into the 
halo. Alternatively, galaxy halos can form instability- 
triggered clouds, resulting in cool g as that is neither 
freshly accreted nor outflowing (e.g. Mailer & Bullock! 
l200i IKaufmann etaDl20"0ll IKeres fc Hernauiso[2009lR . 

This picture of gaseous halo formation has direct 
consequences on observations of halo absorption sys- 
tems. We predict that galaxies above the transition 
mass have drastically smaller covering fractions within 
a fixed co-moving radius, with respect to accreted cool 
gas. While we have presented a small range in val- 
ues for M s h, different models vary in detail. Typically, 
~ 1O U ' 5 ~ 12 ' 5 M0, corresponding to central galaxies 
with stellar masses M star ~ 10 10 - 2 ~ 10 - 8 Af Q at z = 0, 
or M st ar ~ 10 9 8 ~ 10 ; 8 M W at z = 1 (based on abun- 
dance matching, e.g.. lConrov fc Wechsler|[2009l ) . If mod- 
erately high covering fractions are found to exist for 
more massive galaxies (which are not undergoing gas- 
rich mergers), that gas is likely the result of outflows, and 
should show distinct differences in metallicity and kine- 
matics when compared to freshly accreted cool gas (for 
an observable kinemati c signature of co ol accreted gas, 
see IStewart et al.|[20ilT) . We note that iTinker fc Chen! 
(|2008| ) found observational evidence from correlation 

11 To emphasize that a precise value of M s h is uncertain, the 
vertical hashed bar in each panel shows a small range of values, 
M sh = 0.9 - 1 X 10 12 M Q . 

12 While our resolution does not allow us to directly simulate 
these clouds, we do not expect high covering fractions due to 
cloud fragmentation, for standard hot gas profiles ( Kaufmann et al. 
12003 ). 
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Fig. 2. — Covering fraction of neutral hydrogen as a function of the halo virial mass. The left and right panels show Halo 1 and Halo 2. 
Different lines correspond to outer radii of R < 30, 50, 100 co- moving kpc (dot-dashed blue, solid black, dotted red, respectively). Note the 
drastic transition in covering fraction at the cold/hot mode transition of M v - lr 
z ~ 0.8 for Halo 2. 

functions that supports this model of a cold/hot mode 
transition influencing the presence of halo absorption sys- 
tems. 



^ah ~ 10 A/q, which occurs at z ~ 1.3 for Halo 1 and 

transition this sharp should be directly observable 
via metal line absorption system studies. 



5. CONCLUSION 

We have used two high-resolution, cosmological SPH 
simulations as a tool for studying the cool gaseous ha- 
los of galaxies. By creating mock observation sight 
lines, representing cool gas detection via absorption, we 
present the covering fraction of neutral hydrogen in our 
galaxies over time. Our primary results are summarized 
as follows. 



To first order, the covering fraction of accreted cool 
gas is relatively stable at f c ~ 30 — 40% (R < 
50 co- moving kpc, Nm > 10 16 cm -2 ), as long as 
the galaxy continues to accrete cool gas from the 
cosmic web. 



As soon as our simulated galaxies cross the thresh- 
old between cold mode accretion and hot mode ac- 
cretion (when the halo is massive enough to sup- 
port stable shocks at a large fraction of the virial 
radius, M v - n ~ 10 12 M Q ) the lack of cool accreted 
gas results in a suppression of cool halo gas. Within 
~ 500 Myr of reaching this threshold mass, the cov- 
ering fraction drops from 30 — 50% to 5 — 10%. A 



We have used a feedback model without cool gas out- 
flows here, focusing on galaxy halo properties that are 
a natural consequence of cosmological gas accretion in 
LCDM. However, observations have shown that galaxy 
outflows are an abundant phenomenon, and likely play an 
important role in shaping the properties of cool gaseous 
halos around galaxies. We believe future work in com- 
paring observations to a variety of simulations with dif- 
ferent feedback models would prove a valuable tool in 
testing theoretical models of galaxy formation, as well 
as understanding the underlying nature of galaxy halo 
observations. 
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